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ABSTRACT 

We investigate in detail the parameter space of active-sterile neutrino 
oscillations that amplifies neutrino chemical potentials at the epoch of Big 
Bang Nucleosynthesis. We calculate the magnitude of the amplification and 
show evidences of chaos in the amplification process. We also discuss the 
implications of the neutrino chemical potential amplification in the Big Bang 
Nucleosynthesis. It is shown that with a ~ 1 eV i> e , the amplification of its 
chemical potential by active-sterile neutrino oscillations can lower the effective 
number of neutrino species at Big Bang Nucleosynthesis to significantly below 
3. 



PACS: 14.60.Pq, 14.60.St, 26.35. +c, 95.35. +d 
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I. INTRODUCTION 



Neutrino oscillations have been suggested to explain several experimental results and 
if proven true, they will represent a significant step toward physics beyond the standard 
model of particle physics Q. Mixings between active neutrinos {y e , or v T ) and sterile 
neutrinos (hypothetical neutrinos that do not interact with known particles via the strong, 
weak or electromagnetic interactions) are one possible source of neutrino oscillations. The 
most stringent constraints on the parameters of active-sterile neutrino mixings come from 
cosmological and astrophysical considerations In particular, the active-sterile neutrino 

oscillations at the epoch of Big Bang Nucleosynthesis (BBN) have been investigated exten- 
sively and tight constraints have been obtained based on the primordial 4 He abundance 
in our universe. Interestingly, it was recently pointed out by Foot, Thomson and Volkas 
that some parameter space of the active-sterile neutrino mixings can amplify neutrino 
asymmetries (neutrino chemical potentials) so that the previous constraints on active-sterile 
neutrino mixings based on BBN can be alleviated |7]||. Also an asymmetry in the electron 
neutrino sector as large as ~ 0.1 at the time of BBN can change significantly the BBN 
prediction of the primordial 4 He abundance. 

In section 2 of this paper we expand the original investigation of Foot, Thomson and 
Volkas [0, by calculating the parameter space that amplifies neutrino chemical potentials 
and the magnitude of the amplification. But instead of relying on a simplified equation that 
only applies outside the resonant regime and when neutrinos are incoherent, we analyze the 
problem based on the original equations in the density matrix formalism, both analytically 
and numerically. Our analyses reveal many interesting features of the amplification process 
that cannot be revealed by the simplified approach. For example, the neutrino asymmetry 
can be oscillatory long after the initial resonant crossing. There are evidences which suggest 
that the oscillatory asymmetry is chaotic. As a result, although the order of magnitude of the 
final neutrino chemical potential is readily predictable, the sign of the chemical potential is 
very sensitive to the mixing parameters and the input parameters of numerical calculations. 
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This oscillatory behavior and evidences of a chaotic amplification are probed in section 2. 

In section 3, we discuss two implications of our results in section 2 and show how active 
neutrinos as dark matter candidates-having ~ leV mass- can lower the effective number of 
neutrino species in BBN to significantly below 3 by mixing with a lighter sterile neutrino. 



II. FORMALISM AND CALCULATIONS 

Through out the paper, we adopt a unit in which h — c — k — 1. We also use a 
convention to denote the number density of a particle i by N iy and the number density 
relative to its equilibrium value (2((3)T 3 /n 2 for photons, 3/4 of that for electron-positrons, 
and 3/8 of that for neutrinos) by n^. The Hubble expansion rate is H = 7T 2 /M P where 
M p is the Planck mass and T the temperature of the universe. Finally, T$ denotes T in the 
unit of MeV. Since we are only concerned with the era of BBN, we limit our discussion to 
1 < T 6 < 100. 

Mixtures of an active neutrino v a and a sterile neutrino v s can be described by a density 
matrix ||3|| 



Pv 



( \ 

Paa Pas 



P " I + P - S , (i) 



y Psa Pss j 

where a are the Pauli matrices. The number densities of v a and v s in the mixture, relative 
to their equilibrium values, are respectively 

PO + Pz PO — Pz y„s 

n » a = 2 ' Uus = 2 ' 

The evolution of the total relative number density of the neutrino mixture at the epoch of 
BBN is |H 

P = (r^a^a-* ii)){nini-n Ua n Pa ), (3) 

where (r) are reaction rates averaged over a thermal spectrum. Values of (T) are listed in 
table 1 of ref. or ref. f|. The evolution of P is |§ 



V x P + P z - DP_ 



(4) 



where V represents the frequency and the axis of the oscillation in the P-space, and P^ = 
P x x + P y y. The D-term represents the damping of Pj_ due to neutrino interactions which 
constantly reduce a mixed neutrino state into an eigenstate of either v a or v s . 
At the epoch of BBN, 

5m 2 



filTI 2 

V x = —sm2B, V v = 0, V z 



2E 



cos2fl + Kf + K 



(5) 



where 5m 2 and 9 are the usual vacuum mixing parameters, and E is the energy of the 
neutrinos. V£ is the contribution of the matter effect from asymmetries in the background 
plasma |§: 



= V2G F N 7 {L + 0.375 
« 0.13G F T 3 \8L /3 + 2(n 



2(n„ 



\n V a n u ) 



(6) 



where Lq represents the contributions from the baryonic asymmetry as well as the asymmetry 
in electron-positions, and is ~ 10~ 9 . iV 7 is the photon number density. The n u — n D terms 
represent the asymmetries in active neutrinos and thus their non-zero chemical potentials. 
If ^-the chemical potential of v divided by kT-is much smaller than 1, n v — rip « 1.8C, U . 
V£ is the contribution of the matter effect due to a finite temperature |§: 



Vl = -V2G F N 7 [l2MET(n Ua + n Pa )/4M 2 + 12MET/M ] 
= -V2G F NJl2MET{n Va +n Da )/AM 2 



a 



a — ii, t. 



(7) 



It has been shown that eq. (0)-([7D give a good description of neutrino oscillations in BBN 
if the average neutrino energy E ~ 3.151T is inserted in V and if D is thermally averaged 
Pj. Therefore, numerically, 



V 1 



-250GIT 5 , 



-70G|T 5 , 



a 



a: 



e; 



(8) 
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The damping coefficient D, consisting of contributions from both elastic scatterings and 
inelastic scatterings of u a , is |3|,^] 

D fa (1.3 + 0An 1/a + 0.5np Q )GpT 5 , a = e; (9) 
w (0.8 + 0An Ua + 0.5n p JG 2 F T 5 , a = fi,r. 

The reason we leave out n Va and np a explicitly without approximating them to 1 is for 
the convenience of calculating the difference in the coefficient between the v a -v s and i> a -i> s 
oscillations. Since we often compare D to the Hubble expansion rate H, we note 

D « 0.5T 6 3 if. (10) 

The initial condition for eq. ([3]) and (Q) is usually chosen to he P$ = P z = 1 and 
P x = P y = 0, at T init when dominates over 8m 2 /2E, i.e., 

5m 2 cos 26 



T init > 15 



eV 2 



1/6 

(11) 



(for the moment we assume any neutrino asymmetry is negligible). That is, the neutrino 
ensemble consists purely of v w which is a good approximation because V z ^> V x so that V 
is almost aligned with the z-axis. 

As the universe expands and its temperature drops, \V X /V Z \ becomes larger, the ampli- 
tude of the oscillation consequently increases. Eventually, if there is no amplification of 
neutrino asymmetries, V a becomes negligible, and V settles down into its vacuum value. 
During the process, if the mixing has Sm 2 < (v a heavier than v s ), a resonance can occur 
when V crosses the x-axis at a temperature 



T res w 13(16) 



Sm 2 cos 28 



1/6 



MeV for a = e (//,t). (12) 



leV 2 

P and V before and after the resonance are illustrated in figure 1. 

During the oscillation, the interactions between u a and the background plasma play 
two roles. First, the interactions (including both elastic and inelastic ones, represented by 
the D-term) reduce mixed neutrino states into either v a or u s , which effectively damp the 



amplitude of Pj_ and at the same time randomize the phase of the neutrino oscillation. When 
the regenerated v a (and v s but mostly u a ) oscillate into v s (and v a ) again, the portion of v a 
in excess of u s , P z , decreases toward 0. Secondly, the inelastic process-^ a z/ a pair productions 
(the P^z term)-constantly replenishes the number of v a that is being depleted by oscillation, 
maintaining its population as a full relativistic species as long as such pair productions are 
potent (T > 3 MeV for v e and > 5 MeV for or u T ). 

Eq. ©-(0) can be equally applied to the anti-neutrino sector, with notations for par- 
ticles and anti-particles switched and L replaced by — Lq. Apparently, since v a and v a can 
only be produced as pairs, Pq = Pq. 

The relevance to BBN comes at T ~ 1 MeV when the neutron to proton ratio freezes out. 
If a significant population of u s is produced, or a significant asymmetry in v e v e is generated 
through the active-sterile neutrino oscillations, the neutron to proton ratio can be affected 
and the resultant 4 He primordial abundance altered from the standard BBN predictions. 
When neutrino asymmetries are negligible, to be consistent with the observed primordial 
4 He abundance requires 

5m 2 sin 4 29 < 10~ 9 eV 2 (13) 

(The bound on the v e -v s mixing is tighter on the low Sm 2 end. See refs. |§ and H for precise 
constraints.) 

Under conditions that 

|V| >£>> |V|/|V|, (14) 

i.e., the damping of P and the change in V are negligible within one cycle of oscillation of 
P, eq. (Q) can be simplified to the lowest order to 

P X = V X P Z /V Z , P y = 0, P Z = -DV 2 P Z /(V 2 + V 2 ) + P . (15) 

In the epoch of our concern, and in the absence of an amplification of neutrino asymmetries, 
eq. (|i~4l) is satisfied except near the resonance region where \V Z \ ~ \V X \. This is because 
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V? = £>(10 2 )P> and \V a /V a \ ~ H < D at T 6 > 1. We will discuss the case of amplified 
asymmetries later in the section. 

Similarly for anti-neutrinos, the approximate equations are 

P X = V X P Z /V Z , P y = 0, P Z = -DV 2 P Z /(V 2 + V 2 )+P (16) 

under conditions that |V| ^> D ^> |V|/|V|. Assuming there is no asymmetries in neutrinos 
other than the oscillating v a v a sector, 

V Z = V Z - 2f3(AP z + 8L /3) =V - (3(AP Z + 8L /3) (17) 

where (3 = 0.375y/2G F N y » 0.13G F T 3 , V = -5m 2 cos 29 /2E + V? and AP 2 = P z - P z = 
2{n Va - n Da ). It is also noted that D - £> « O.OSAP^T 5 < P>AP 2 . 

The asymmetry in the ^Pq, sector can then be described by AP Z which to its lowest 
order satisfies 

2 T/2 



Ah = ^(^5 - ^) " ^T7^AP 2 + (P - D)J^R 



~v x 2 + v? v 2 + v 21 v? + v? v 'v* + v* 

DV 2 \ 4V /3(AP 2 + 8L /3)P 2 



1 }AP Z . (18) 



+ [V - I3(AP Z + 8L /3)] 2 { {V* + [V + P(AP Z + 8L /3)] 2 }AP 2 

The equation resembles equation (15) of Foot, Thomson and Volkas pj except that their 
equation omitted the second term and L , and assumed P z to be 1. 

When Vq < 0, eq. (|18D is a damping equation for |AP Z | > 8L /3, and no amplification of 
|AP 2 | to 3> 10~ 9 can occur. This rules out i/ a -i/ s mixings with 5m 2 > (i/ s heavier than u a ). 
Only when 5m 2 < and Vq switchs to a positive value (resonance crossing) before T ~ 1 
MeV does an amplification of AP 2 become plausible. 

For mixings with 5m 2 < 0, Vq still starts out negative at high temperatures 3> T res . Any 
\AP Z \ > 8L /3 will be damped toward an asymptotic value such that AP 2 — ► 0. Thus 

8V 2 L /3 



AP Z -> --L + 



4Vq/3P, 



-3^0 (19) 



when I Vnl 3> IVU. 



When the system enters the resonant regime, eq. ( I8[) does not apply. We have to go 



back to the original equation (Q) and its anti- neutrino counterpart, which give 
AP X = -V AP y - (3(AP Z + 8L /3)(P y + P y ) - PAP,, 

AP y = V AP X + (3{AP Z + 8L /3)(P X + P x ) - V X AP Z - DAP y , (20) 
AP Z = V x AP y . 

where AP X = P x — P x and AP y = P y — P y . If we don't want v s to be brought into equilibrium, 
we have to restrict our discussions to the parameter space that satisfies eq. (|13|). Then: (1) 
the resonance crossing is non-adiabatic; (2) D has to be small enough compared to the 
time scale of the resonance crossing so that most of the neutrinos do not scatter during the 
crossing. The non-adiabatic change in V leads to a coherent oscillation of P around the new 
V (see figure 1(b)). The impotency of scatterings enables this coherency to be maintained 
throughout the resonant regime and beyond until D > |V|/|V| so that the interactions have 
enough time to randomize the phases of neutrinos again. It also implies that the D-term in 



eq. fl20D can be dropped in resonance. 

When the temperature was high above T res so that eq. fllTf ) applies, AP X ps V x AP z /V z , 
AP y and AP y are approximately 0. But in and right after resonance, the approximation 
breaks down because of the rapid change of V z and V z . Instead AP y becomes oscillatory 
with a frequency of ~ V z (figure 1(b)). \AP Z + 8L /3| will quickly be of order 10~ 9 , but 
a more important question is whether \AP Z + 8L /3| (and thus |AP Z |) can be amplified to 
^> 1CT 9 . Assuming that \(3AP Z \ <^ Vo, the amplitude of the oscillating AP y will be of order 
~ \P Z V X {V Z X — Kr 1 )! ~ \2V x /3AP z /Vq \ (from now on P z ~ 1 and L is dropped for simplicity 
because \AP Z + 8L /3| ~ |AP 2 |). The amplification of \AP Z \ depends on whether AP 2 has 
enough time to change AP 2 by a factor of more than 1, i.e., 

win 



> V . (21) 



Since Vo is a changing quantity, the condition of amplification depends on which Vq to choose. 
For a crude estimate, a reasonable choice is the Vq at the time when AP y oscillates one cycle 
since the resonance (so that eq fl2l|) is meaningful). So 
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V o-Vo-V -^^.\^i. (22) 
Solving the equation assuming E = 3.151T yields Vq ~ 10~ 2 |5m 2 | _a25 ■ \5m 2 /2E\. Thus the 



condition of amplifying AP 2 , i.e., eq. fl21|) , is 



|5m 2 r 5 / 12 sin 2 2#> 1(T 12 . (23) 

The growth of |AP 2 | is limited once \@AP Z \ > V (when |AP 2 | > |5m 2 /9eV 2 |T 6 " 4 ). 
Because at this moment, the amplitude of AP y becomes ~ \2V X /(3AP Z \, and the amplitude 
of AP 2 /AP 2 becomes ~ 2V; 2 //3AP 2 2 , proportional to AP7 2 . 

An interesting feature of AP 2 , shown both from eq. ([H]) and numerical calculations, 
is that it keeps oscillating below the resonant temperature (figure 2). This is a direct 
consequence of the coherent oscillation of P and P. As a result of the oscillation, the sign 
of AP 2 flips (so does V z and V z ) unless the change of AP 2 in each cycle is smaller than the 
amplitude of AP 2 itself, i.e., 

|AP -' > wk ' lpAPrK (24) 

Since for parameters that satisfy eq. (|23|), the amplitude of AP Z > (|5m 2 |/9eV 2 )T 6 ~ 4 , a rough 
realization of the condition that AP 2 will not be oscillatory between positive and negative 
values is 

IS a (*) 

At T 6 = 1, AP Z will not be oscillatory if \5m 2 \/ sin 2 29 > 9eV 2 . 

Once AP 2 stops flipping its sign, |V|/|V| and |V|/|V| decrease dramatically so that 
the damping P-term becomes more and more important. Eventually eq. fll8| ) reapplies. 
Assuming P z ~ 1 and neglecting L and V x in the denominator eq. (|I8|) yields 

W (3 



V 

AP w D 



(V -(3AP z y 



(V + (3AP Z ) 2 



AP Z . (26) 



Since Vq > after the initial resonance, this is an amplification equation if 



4VW3 , .. , „., 



(V + (3AP z )z 
The amplification will not stop until 



> 1, or 5m 2 < 10 2 eV 2 . (27) 



p V x WoP < „ (2R) 

(V + {3AP z y(V -f3AP z y~ ' 1 ' 



which occurs at \/3AP z \ > V ~ |<5m 2 /2P| and 

/l8T 7 sin 2 2#\ (5m 2 /2E\ 4 



< I- (29) 



V |5m 2 | / V /3AP 2 

At Tg ~ 1, nearly all the mixing parameters that show no oscillatory AP 2 have 
18T 7 sin 2 2fl/|<5m 2 | < 1, so \(3AP Z \ is limited to ~ \5m 2 /2E\ and 
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at a temperature of ~ 1 MeV. This limit is confirmed by our numerical calculations and 
is similar to that of Foot, Thomson and Volkas @ based on their simplified equation. For 
\Sm 2 \ > 10eV 2 , however, since AP Z has to be much smaller than P z ~ 1, AP Z ~ 0.1. 

The oscillatory behavior of AP Z after resonance is illustrated in figure 2 (a), (b) and (c), 
for three different v e -v s mixing parameters. Each graph is the result of several millions steps 
of integrations of eq. (|) and eq. (|J) by adaptive Runge-Kutta method, with an error of less 
than 10~ 10 in each step. In figure 2 (a), the mixing parameters do not satisfy equation ( pop 
at T 6 = 1, so AP Z is still oscillatory at 1 MeV. In figure 2 (b) and (c), |AP Z | settles down to 
~ |<5m 2 /9eV 2 |T 6 -4 at T 6 > 1, in line with our estimate eq. (|30|). Our numerical calculations 
also show that although the final settle-down value of |AP Z | is predictable, the sign of AP Z 
seems random among different parameter choices. For example, in case (b), 5m 2 = — 10 -2 
eV 2 and sin 2 28 = 10~ 4 25 , a small change of sin 2 29 to l(J -4 - 1875 yields an opposite sign of 
the final AP Z . The sign of AP Z can also be flipped by slight changes in the initial L (as 
tiny as 0.01%) and calculational parameters, such as a different error control (from 10~ 10 
to 5xl0~ 10 in our example), or a step size, or even a slightly different relation between 
the average neutrino energy and the temperature (from E = 3.151T to E — 3.150T in 
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our example). This is due to the large number of oscillations of AP Z before it approaches 
one of the two possible values, so that the final AP Z is very sensitive to the input and 
calculational parameters. Such behavior is not so significant in the case of 5m 2 = — 1 eV 2 
and sin 2 29 = 10~ 8 , because the number of oscillations of AP Z is small (figure 2 (c)). In this 
case, using E = 3.150T the evolution of P and P traces the evolution of P and P using 
E = 3.151T very well (in a sense that their difference is obviously still a perturbation at 
lower temperatures). This is also true if we change the initial L by 0.01% (although the 
resultant perturbation can be several percent in the oscillatory epoch of AP Z ). Nevertheless, 
a small change of sin 2 29 to io~ 8 - 0625 still flips the sign of the final AP Z . 

Figure 2 (a) and (b) suggest a chaotic behavior in the epoch of oscillating AP Z . Figure 
2 (c) might be intrinsically chaotic too but the time scale of the oscillatory epoch may 
be too short for such behavior to show up. To find more evidences of chaos, we try to 



determine the Lyapunov exponents of the system fllQ"| . Assuming 0(t) = (AP X , AP y , AP Z ) 
represents a solution of eq (|20"D, we investigate the behavior of a nearby solution (j)(t)+5(j)(t) = 
(AP X , AP y , AP Z ) + (SP X , 5P y , SP Z ) {5P X and SP y are assumed to arise from P x and P y only 
so that we can carry out our analysis). The evolution of the 5(p(t) is obtained by linearizing 
eq. 



50 = M5(j) 



(31) 



where matrix 



M 



V 



~D -V - f3AP z -P(Py + Py) 

-V -(3AP Z -D - Vx + f3{Py + Py) 
OK 



D. 



(32) 



(33) 



The eigenvalue of M, called the Lyapunov exponent, is 

(3{Py + Py){V +(3AP z ) 

(3(P X + P x ) - V x 

If A is positive, nearby solutions will depart exponentially in phase space within a timescale 
of A -1 . In our problem, since P x + P x , P y + P y and AP Z are all oscillatory, so is A. A crude 
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analysis is to plug in the amplitudes of P x + P x , P y + P y and AP Z , to see whether A can 
be sometimes positive. In the limiting case of \/3AP z \ <C V , P x + P x ~ P y + P y ~ V^/Vo, 
the resultant A is ~ (3Vq/{(3 — Vb) — -D which can certainly be positive if Vo > and 
1 5m 2 1 <C 9TgeV 2 , conditions satisfied after the initial resonance. In the limiting case of 
\(3AP Z \ > V , P x + P x ~ Py + P y ~ V x /f3AP z , the resultant A is ~ PAPJ(1 - AP Z ) - D, 
which can again be positive (remember AP Z <C 1 and (3AP Z > > f). Therefore, 
our crude analysis shows that A can at least be positive within a timescale of order A -1 
intermittently!] In other words, the system is not a classical textbook example of a chaotic 
system. 

We are ultimately interested in the region of mixing parameters that amplifies neutrino 
chemical potentials, and the size of the amplification. In figure 3 (a) and (b), we plot the 
parameter space allowed by BBN that amplifies neutrino chemical potentials. The boundary 
to the right which excludes parameters that bring v s into equilibrium is adopted from Shi, 
Schramm and Fields || . The boundary to the left that distinguishes parameters that amplify 
neutrino chemical potentials from those that do not is based on our numerical calculations 
(smoothed), and agrees with our analytical estimate eq. ( |2"3"D within an order of magnitude. 
The lower cut on \Sm 2 \ is determined by requiring T res > 1 MeV. The upper cut on \Sm 2 \ is 
dictated by laboratory bounds on u e mass in the v e -v s mixing case, and by eq. ( P7f) as well 
as cosmological considerations |TT| in the Ufj,(u T )-u a mixing case. The boundary that singles 
out parameters that have oscillatory AP Z at 1 MeV is plotted according to eq. (|25|) which 
is confirmed by our numerical calculations. 

We note that our numerical calculation of 5m 2 = — leV 2 , sin 2 26 = 10~ 8 v e -v s mixing 



In the analysis, we have chosen a particular set of nearby solutions, namely those having devia- 
tions in P x and P y but not in P x and P y . The other extreme choices, in which 5P X or SP y arises 
only in one particle population but not in its anti-particles, merely change the sign of /3AP Z in the 
nominator of eq. ([£]), thus do not affect the general behavior of A. 
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yields an opposite sign of AP Z from that of Foot, Thomson and Volkas 0. But this may not 
be surprising due to the chaotic feature of the system, that different signs of AP Z may arise 
from different initial Lq, or even different choices of integrators, different errors or step sizes. 
It is also noted that the simplified equation in ref. |7j (corresponding to eq. flI8|) without the 
second term and Lq) is not suited for investigating the behavior of AP Z in the resonant regime 
and in the epoch of oscillatory AP Z thereafter. Finally, we note that the first calculation of 
the neutrino asymmetry done by Enqvist et al. [|12|] shows an oscillatory asymmetry down 
to Tq ~ 1, because their parameter choice, Am 2 = — 10 -5 eV 2 and sin 26 = 10"2, does not 
satisfy eq. (f25[). 



III. IMPLICATIONS 

We concentrate on two implications of a neutrino asymmetry as large as eq. ( |30| ) in BBN. 

The first is on other active-sterile neutrino oscillations in BBN [|]|§. When neglecting 
neutrino asymmetries, large areas of parameter space of active-sterile neutrino oscillation are 
ruled out based on the argument that the sterile neutrino cannot be significantly populated 
so as to violate the primordial 4 He abundance observation. The forbidden areas include the 
large angle v^-Vg mixing with 6m 2 ~ 1CT 2 eV 2 which can solve the atmospheric neutrino 
problem, and the large angle v e -v s mixing with 5m 2 ~ 10~ 5 eV 2 which can solve the solar 
neutrino problem. This argument, however, no longer stands when a neutrino asymmetry 
as large as in eq. (^) is in place. For example, if v T mixes with a lighter v s with \Sm 2 \ ^> 
l(T 2 eV 2 and an angle in the shaded region of figure 3 (b), the v T v T asymmetry amplified 
by the v T -v s oscillation can be large enough to suppress the v^-v,, oscillation from the v^-v,, 
mixing solution to the atmospheric problem ||. Similarly, if v T or z/ M mixes with a lighter 
v s with 1 5m 2 | ^> 10~ 5 eV 2 and an angle in the shaded region of figure 3 (b), the amplified 
neutrino asymmetry can be large enough to suppress the v e -v s oscillation originating from 
the large angle v e -v s mixing solution to the solar neutrino problem. The suppression is in 
place even in the epoch of oscillating AP Z , because although an oscillating AP Z constantly 



13 



drives the other active-sterile neutrino oscillations through resonance, the time of resonance 
crossing is too short to allow any significant oscillation. Thus both these two solutions to 
the atmospheric neutrino problem and the solar neutrino problem ruled out previously may 
still be viable if an active neutrino (more massive than or v e respectively) mixes with a 
lighter sterile neutrino with parameters in the shaded region of figure 3 (b). 

The second implication is on the primordial 4 He abundance itself. Besides the number 
of neutrino species, the primordial 4 He abundance is also affected by a non-zero chemical 
potential in the v e v e sector at T ~ 1 MeV. The mechanism is that the asymmetry in 
v e v e changes the neutron/proton conversion rates, thereby changes the freeze-out time of 
the neutron to proton ratio as well as the ratio itself. A £ Ve (the v e chemical potential 
divided by kT) of order 0.1 can induce an appreciable change in the prediction of the 
primordial 4 He abundance 0. When £„ e < 1, Y « Y(^ e = 0) - 0.234f„ e and 
AP Z = 2{n Ve - n Pe ) « 3.6£„ e . So 

Y = Y - 0.065AP,. (34) 

The comparison with equation Y = Y(N U = 3) + 0.012(iV I/ — 3) (where N v is the effective 
number of neutrino species in BBN) [[T5[] indicates that AP Z ~ 0.1 in the z/ e z/ e sector cor- 
responds to roughly —0.55 neutrino species, and therefore has a significant impact on the 
predicted 4 He abundance. 

There are two ways to generate a £„ e of order ±0.1 by active-sterile neutrino oscillations. 
The direct way is to have a ~ 1 eV v e mix with a lighter v s (figure 3 (a)). If the atmo- 
spheric neutrino problem and the solar neutrino problem are to be solved by active neutrino 
oscillations, this implies that all three active neutrinos are almost degenerate with a mass 
of order 1 eV. This will be consistent with supernovae nucleosynthesis constraints p!6|Ji7[| 
and compatible with the controversial LSND result [fT8] , [3~9[1 if the claimed detection of v^-Ve 
oscillation solves the atmospheric neutrino problem ||20|| . Laboratory experiments limit the 
mass of v e to less than 5 eV [BTJ. If v e is a majorana neutrino, its mass is further limited to 



less than about 1 eV 22 
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The indirect way of generating a significant £„ e is to have u T (or v^) mix with a lighter 
v s with 10 2 eV 2 > \5m 2 \ > 1 eV 2 and a desired angle, and transfer the asymmetry in the 
v T v T {PyVfj) sector into v e v e by a v e -v T (z/ M ) mixing. But to yield an asymmetry of order 0.1 
as well in v e v e ^ the transfer has to be efficient. Take the v e -v T oscillation as an example, the 
mixing has to satisfy 

D '(y)j *H, (35) 

where D', V' x and V' z denote the counterparts of D, V x and V z in the v T -v s oscillation. 
Approximately 

D' ~ 5G 2 F T 5 w T 6 3 #, 
x 6.3T 

V'z « cos2 ^ " 180G|T 5 + 0.13G F T 3 (n, e - n p J - 0.l3G F T 3 (n„ T - n Pr ), (36) 

where 5M 2 and 0' are the vacuum mixing parameters of the v e -v T mixing. An efficient 
transfer of asymmetry means that (n Ue —n De ) ~ {n Vr —n PT ) ~ |<5m 2 |/18Tg. If the mass of v e 
is much lighter than v T) \5M 2 \ pa \5m 2 \. Then (D/H)(V X /V Z ) 2 ~ T 6 3 sin 2 26>' at temperatures 
approaching 1 MeV. So if i/ r has a cosmologically interesting mass, Am 2 ~ 20-100 eV 2 , 
| AP Z \ pa 0.1 in v T v T can be reached at T pa 2-3 MeV by the v T -v s mixing according to eq. fl30"l) , 
which can efficiently transfer into an asymmetry of similar order in v e v e if v T mixes with v e 
with sin 2 26' ~ 0.1 and |5M 2 | ~ |<5m 2 |. This required mixing between u e and u T lies near the 
edge of current lab limits on the v e -v x mixing f23|| , and may be testable in the near future. 



Based on supernovae nucleosynthesis arguments, however, the required v e -v T mixing is ruled 
out (although supernova models are uncertain to some extent) |16| , [17| . The above analyses 
applies similarly to the v^-Vg and v^-v e mixings, but the required mixing between v e and 
to transfer asymmetries efficiently has already been ruled out by laboratory experiments 
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Of course, the indirect way of transferring asymmetries from v T v T (or v^v^) into v e v e 
works if they have almost degenerate masses, ~ leV. In this case the required Ve-v^ (u T ) 

15 



mixing will not be ruled out by laboratory experiments or astrophysical considerations. 

Recently, Hata et al. |2|] suggested a possible "crisis" in BBN because the predicted 4 He 
abundance from the standard BBN, coupled with predictions of 3 He and D from generic 
chemical evolution models, is too high to be consistent with the observed abundance, unless 
the number of effective neutrino species N v at the time of BBN is less than 2.6. This is at 
odds with popular beliefs and most theoretical models which assume all u e , and v T to be 
lighter than ~ 1 MeV and therefore N v > 3. But as seen above, if an asymmetry of order 
AP Z ~ 0.1 arises in the v e v e sector from active-sterile neutrino oscillations, this potential 
"crisis" can be solved. Of course, the mixing parameters of the active-sterile oscillations 
have to be right to yield a positive AP Z instead of a negative one. An interesting note is 
that the neutrino mass required to solve this "crisis" is ~ 1 eV, which qualifies neutrinos 
as dark matter candidates. A v e mass of ~ 1 eV is within a factor of 5 below the current 
limit on v e mass, and right on the edge of detection limit if v e are majorana neutrinos. 
If not introducing more sterile neutrinos, the solar neutrino problem and the atmospheric 
neutrinos have to be solved by mixings among the three active neutrinos, therefore requiring 
their masses to be almost degenerate. A model of three neutrinos almost degenerate in mass 
is favorable in forming structures in the universe [2S|, but may not be the most natural 
theoretical model so far. 



IV. SUMMARY 

In summary, we have calculated the parameter space of active-sterile neutrino mixings 
that amplifies neutrino chemical potentials, and the size of the amplification. Results are 
summarized in figure 3. By exploring the sensitivity of the amplification to the initial 
condition, the mixing parameters and the calculational parameters, and by analyzing the 
Lyapunov exponent of the system, we showed evidences that the amplification process is 
chaotic. We have also discussed the implications of our results on BBN. It was shown that 
a u e chemical potential of order O.lfcT could be achieved by either a mixing between ~ 1 



16 



eV v e and a lighter sterile neutrino, or a mixing between a ~ 1 eV ^ (or u T ) and a lighter 
sterile neutrino coupled with a mixing between almost degenerate u e and i/ M (or u T ). Such 
a chemical potential in v e can lower the effective number of neutrino species to significantly 
below 3. 
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Figure Captions: 

Figure 1. Illustrations of evolution of P and V for 5m 2 < 0. (a) before resonance; (b) right 
after resonance. 

Figure 2. The solid lines show the evolution of AP Z vs. the temperature of the universe, 
for v e -v s mixing. Asymmetries less than 10~ 9 are ignored. The dash line shows \AP Z \ = 
|5m 2 /9eV 2 |T 6 - 4 . L = 10~ 9 . (a) 5m 2 = -1(T 4 eV 2 , sin 2 26 = 1(T 5 ; (b) 5m 2 = -1(T 2 eV 2 , 
sin 2 26 = 10" 4 - 25 ; (c) 5m 2 = -1 eV 2 , sin 2 26 = 10~ 8 . 

Figure 3. Regions in between the two thick lines represent allowed mixing parameters that 
amplify neutrino chemical potential. The magnitude of AP Z at 1 MeV is shown. The region 
at the lower right noted with U AP Z Osc." has an oscillating AP Z at 1 MeV. For |<5m 2 | > 10 
eV 2 , AP Z is limited to ~ 0.1. (a) v e -v s mixing [y e heavier than i/ s ); (b) u^{u T )-u s mixing {v^ 
or v T heavier than v s ). 
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